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We investigate chaotic, memory and cooling rate effects in the three dimensional Edwards- 
Anderson model by doing thermoremanent (TRM) and AC susceptibility numerical experiments 
and making a detailed comparison with laboratory experiments on spin glasses. In contrast to the 
experiments, the Edwards-Anderson model does not show any trace of re-initialization processes 
in temperature change experiments (TRM or AC). A detailed comparison with AC relaxation ex- 
periments in the presence of DC magnetic field or coupling distribution perturbations reveals that 
the absence of chaotic effects in the Edwards-Anderson model is a consequence of the presence of 
strong cooling rate effects. We discuss possible solutions to this discrepancy, in particular the small- 
ness of the time scales reached in numerical experiments, but we also question the validity of the 
Edwards- Anderson model to reproduce the experimental results. 



I. INTRODUCTION 

One of the most characteristic effects in disordered or glassy systems in their non-stationary regime is the presence of 
aging. The response of the system stiffens with age showing that it depends on all the previous history through e.g. the 
waiting timeum. Experimentally, thisi-pkenomenon is well documented through magnetization relaxation experiments 
and AC susceptibility measurementsoD. Despite of the different experimental procedures needed for both type of 
measurements, magnetization relaxation and AC susceptibility give similar information regarding the aging behavior 
and its waiting time dependence. 

On top of all these non-equilibrium phenomenology, recent dynamical experiments in spin glasses show very peculiar 
chaotic (also called rejuvenation), memory as well as cooling rate effectsaia. These effects are thought to be the 
signature of the spin glass state being much different to those found in usual ferromagnets or other disordered 
systems!! The most unusual experimental result in spin glasses is the absence of cooling rate effects. The approach to 
equilibrium at a given temperature after cooling from high temperatures is not influenced by the whole cooling history 
at higher temperatures but only by the time spent at the last temperature in the thermal history. Experimentally, 
rejuvenation or chaotic effects in spin glasses are measured in a clear way by doing AC measurements. An alternating 
magnetic field of frequency u — -p where P is the period is applied to the sample and both components of the AC 
susceptibility (the in-phase x' and the out-of-phase x") are measured. 

Although the major part of these measurements have been done on insulating spin glasses they are common also to 
metallic spin glasses leading to the question whether these effects are also present in the most well known theoretical 
models. Despite of the large amount of theoretical work devoted to aging effects in glasses and spin glasses there is still 
no convincing and final explanation for the origin of these peculiar chaotic and memory effects. The comprehension 
of these effects will certainly provide a clue to the understanding of the nature of the glassy state. 

Because memory and chaotic (or rejuvenation) effects are intrinsic to spin glasses (metallic and insulating) it is 
important to understand whether models for spin glasses are able to reproduce the experimental results. It is widely 
accepted that the Edwards- Anderson model contains the main features observed in real spin glasses. The purpose 
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of this paper is to present a detailed and critical study of these phenomena in the Edwards- Anderson model in 
three dimensions. This is not a simple matter to address. Despite of the large amount of numerical studies on 
equilibrium and non-equilibrium phenomena there is no clear evidence that the Edwards- Anderson model reproduces 
the main results found in experiments. Note that, e yea t he question whether there is or not phase transition in the 
3D Edwards- Anderson model is still not fully scttlcdHO. 

The purpose of this paper is to present a numerical investigation, fully experimentally oriented, of the non- 
equilibrium behavior of the three dimensional Ising spin glass with special emphasis on recent experiments where 
memory and chaos effects where found. This question is of the out-most importance concerning modeling. If some 
experimentally observed results are missing in any theory then we must understand whiz_. There have been several 
investigations in the literature devoted to this subject, but still a clear answer is missingE3~lij. Here we will provide 
a complementary investigation to results already published, emphasizing the experimental results and comparing 
different types of experiments. In particular, our main effort will be devoted to investigate thermoremanent and AC 
numerical experiments. Although thermoremancijtt|Studies have been largely considered in the past there are very few 
numerical investigations devoted to the AC topict3. 

The paper is divided as follows. In Section |l| we introduce the model as well as the dynamical procedure. Section III 
discusses the two type of measurements we have done: magnetization relaxation and AC numerical experiments. 
Sections IV and [v| present a detailed investigation of memory and chaotic effects with thermoremanent and AC 
experiments respectively. Finally we present a discussion of the results. 



II. THE EDWARDS-ANDERSON MODEL AND SOME DETAILS OF THE SIMULATION 

The Edwards- Anderson modelEl was proposed in the early 70 's as the simplest model which contains the main 
ingredients relevant to explain the spin-glass phenomenology. In particular, it displays a phase transition characterized 
by the onset of freezing in spin-spin correlations and a divergent non-linear susceptibilityEII. The model is defined by 
the following Hamiltonian 

v 

H = — JijCTiUj — h y j <jj , (I) 

(id) 1=1 

where the indices i, j run from 1 to V, the <7j are Ising spins and the pairs identify nearest neighbors in a finite 
dimensional lattice. The exchange couplings Jij are taken from a random distribution. To avoid degeneracy of the 
ground state the simplest choice is a Gaussian distribution with zero average and finite variance, 

The model is defined in any number of finite dimensions although our main concern here is the three-dimensional 
case where there is a spin glass transition at finite temperature T c ~ 95AE3 Hereafter, unless differently specified, 
we will consider A = I without loss of generality. 

Monte Carlo simulations of (|l|) use random updating of the spins with the Metropolis algorithm. A spin is randomly 
chosen and its value changed with the proper probability. Dynamical experiments use very large lattices (typical sizes 
are in the range L = 20 — 100) with negligible finite-size effects for the largest sizes (L = 64 for magnetization 
relaxation experiments and L = 100 for AC experiments). Here we present two classes of different but related 
experiments. Magnetization and correlation relaxation simulations have run on a special purpose machine APEE3 for 
sizes 64 3 and averaging over 10 or 100 samples. AC experiments were run for a single sample on a Linux cluster of 
PC's for sizes L = 64, 100. 

Before presenting the results it is convenient to discuss the fidelity of the EA modcL-to real spin glasses. Clearly, 
the EA model is an idealization of the real microscopic interaction found in spin glassesE-J Spin glasses are commonly 
distinguished into two large classes: metallic and insulating. Metallic spin glasses are diluted magnets where a metallic 
host matrix is doped with some ferromagnetic impurities (for instance AgMn, AuFe, CuMn). In these systems spin 
interactions are due to indirect exchange and mediated through conduction electrons (the RKKY interaction). Metallic 
spin glasses are then diluted magnets where site disorder induces frustrated short-ranged interactions (decaying like 
1/r 3 with r being the distance between impurities). Insulating spin glasses are much different. In this case, exchange 
interactions are usually antiferromagnetic between neighbor spins but dilution and defects lead to a strong frustration. 
Apart from the different microscopic origin of the frustrating interaction, spins are really Hciscnbcrg-likc and the Ising 
behavior arises from the uniaxial anisotropy present in these type of systems. Because anisotropy is usually strong 
and the local rotational symmetry of Heisenberg spins is broken, a treatment taking pure Ising spins turns out to be 
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a good approximation!^. Having in mind these limitations, the Hamiltonian (|l]) is the simplest model which contains 
disorder and frustration, the two ingredients commonly found in real spin glasses. 

III. MAGNETIZATION RELAXATION AND AC EXPERIMENTS 

There are two alternative but equivalent ways to experimentally investigate non-equilibrium phenomena in spin 
glasses: magnetization relaxation experiments and AC measurements. A very complete description of these methods 
can be found inu. Here we only remind the main results. 

• Magnetization relaxation (TRM) experiments 

Relaxation measurements are done applying a uniform magnetic field and measuring the decay of the ther- 
moremanent magnetization (hereafter referred to as TRM), equivalently, the growth of the zero-field cooled 
magnetization. The typical experiment consists in the following. A sample is fastly quenched below the spin 
glass transition temperature for a time t w (i.e. the waiting time). Then a uniform small magnetic field h is 
applied and the growth of the magnetization measured, 

I v 

X(tw,t w +t) = —^ai(t w +t) . (3) 

i— 1 

In the linear response regime Eq.(||) can be written as 

x(tw,t w +t) = J R(t w +t,s)ds , (4) 

where R(t, s) is the response function which gives the change of the magnetization SM at time t when a pulse of 
the magnetic field Sh is applied at previous time s. In spin glasses aging manifests by the fact that x{tw> t w + 1) 
shows a strong dependence on the value of t w . In general, one finds the following decomposition 

X(tw,t w +t)=Xst(t)+Xag(t w ,t w +t) , (5) 

where Xst and Xag are respectively the stationary and aging parts. Experimentally the aging part approximately 
scales with the waiting time t w in the following way 

Xag(tw,t w +t) = f(t/t w ) , (6) 

although systematic deviations from this scaling behavior have been observed. This point will be discussed later 
on. 

Related to magnetization another quantity of interest which can be numerically investigated in simulation are 
two-time correlations. These quantities are difficult to experimentally measure in spin glasses but very easy to 
compute in simulations. They are defined by 

1 v 

c(tw,t w + 1) = — <Ji{t w )(Ji{t w + 1) . (7) 

i— 1 

Again, in the non-equilibrium regime Eq.(Q) can be decomposed in two pieces, a stationary part plus an aging 
part 

C(t w ,t w + t) = C s t(t) + C ag (t w ,t w + t) . (8) 

Similarly to the magnetization, the aging part of the correlation is approximately described by the following 
scaling behavior 

C ag (t w ,t w +t) = g(t/t w ) , (9) 
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again with systematic (but small) deviations respect to it. The stationary part of the correlation and magneti- 
zation are related through the fluctuation-dissipation theorem (FDT) 

Although in mean-field spin glasses a more general relation seems to be valiofi. It links response and correlation 
functions also in the off-equilibrium regime through 



X[C] = -T dx{s ' t] 



dC(s,t) 



(11) 

C{s,t)=C(t w ,t w +t) 



where the fluctuation-dissipation ratio X depends only on the correlation function in the large times limit 
(t,t w — > oo ). In the quasi-equilibrium regime (t < t w ) we have that X — 1 and we recover the usual FDT. 
In the aging regime (t > t w ) the ratio is smaller than one X < 1 and it can be interpreted as a larger 
effective-temperature T c g ~ T/X. In finite-dimensional spin glasses the validitv-of Eq.(|ll|) has been numerically 
checkedEJO and it has been related to the equilibrium distribution of overlaps^!. 

• AC measurements 

In these experiments an oscillating magnetic field h(t) = ho cos(2Trujt) of frequency ui = -p, where P is the 
period, is applied to the system and the magnetization measured as a function of time 

M(f) = M cos(2tt^ + 4>) , (12) 

where Mo is the intensity of the magnetization and 4> is the dephasing between the magnetization and the 
field. The origin of the dephasing is dissipation in the system which prevents the magnetization to follow 
the oscillations of the magnetic field. From the magnetization one can obtain the in-phase and out-of-phase 
susceptibilities defined as 

, _ Af o cos(0) _ 2 f Q P M(t)cos(2irLut)dt 
X h ~ ho ' (13) 

„ _ M sin(^) _ 2/ P M(f)sin(27T^)^ 

X i i v ) 

ho h 

The dephasing <j) measures the rate of dissipation in the system and is given by 

tan(0) = 2L . (15) 
X 



In numerical simulations the in-phase and out-of-phase susceptibilities are computed by averaging the right-hand 



side in Eqs.(13),(|14|) over several periods P = — . This means a very large measurement timeJbr low frequencies 
for both experiments and simulations. In the non-equilibrium regime the AC susceptibilityE3 depends on both 
the waiting time and the frequency. On general grounds one expects that 

X(v,t) =Xai(w) +Xag(v,t) , (16) 

where the aging part of the AC susceptibility approximately satisfies a scaling behavior 

X« fl (w,t) ~h(ut) . (17) 



Both types of measurements give equivalent information about the relaxation dynamics but in different time sectors. 
As discussed in ref.EI, TRM experiments give information on time scales ranging between the two limits t <C t w and 
t 3> t w . For AC experiments the frequency uj corresponds to the inverse of the observation time t b s (note that in 
TRM experiments after switching the field we have that t b s = t) while the age t a corresponds to the total elapsed 
time t + t w . In AC experiments in order to get reliable results on \' and x" one needs to average over several periods 
of the field, while keeping the age of the system more or less unaltered (otherwise the two limiting regimes would 
mixed and the results would be unclear). This is possible only if Lut a 3> 1, which imply t <ti t w . Consequently, is, 
AC measurements one is able only to explore the beginning of the aging regime, also called quasi-stationary regimeES 
which represents a smaller time window than in TRM experiments. 
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IV. MEMORY AND CHAOS IN CORRELATION AND RESPONSE FUNCTIONS 



In this Section we perform a study of memory and chaos effect in spin glasses measuring correlation and response 
functions. We always take the measurements from 10 samples of a 64 3 system (unless differently specified). We closely 
follow the experimental procedure on what concerns temperature changes and we keep the ratio between time scales 
entering in the simulation similar to the experimental ones (with the same limitations of sizes, magnetic fields and 
absolute time scales as already discussed). 



A. "Cooling and stop" experiment 
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FIG. 1. The "cooling and stop" experiment in the EA model gives no evidence for such strong chaos and memory effects 
experimentally found on real samples. Here we use 100 samples of a 64 3 system, cooling rates inversely proportional to P and 
probing time scales O(P). In the inset we show the reference curves measured without any stop during the cooling. 



The experiment is performed in the following way (for more details the reader is address to the original paperO). 
Starting from the high temperature phase, a spin glass sample is cooled at a fixed cooling rate into the glassy phase. 
When a temperature T* ~ 0.8T C is reached, the cooling process is stopped and the sample is let to relax for a long 
time. This relaxation produces a decrease in the susceptibility (in both x' and x") with respect to the reference curve 
(obtained with the same constant cooling rate and without any stop). After that long time the cooling process is 
continued down to a low temperature and finally the sample is heated back again at a constant heating rate (equal 
to the cooling one) and without any stop. 

There are two relevant results in this experiment. First, as soon as the cooling process is started again after the 
stop, the susceptibility merges rapidly with the reference curve, quickly "forgetting" the thermalization work done 
near the temperature T* (chaos effect). Second, when the sample is heated back at a fixed heating rate and without 
any stop, the susceptibility closely follows the cooling curve and it goes through the dip at T* (memory effect). 

In our simulations we do an analogous experiment, measuring correlation functions instead of susceptibilities. We 
divide the time of the experiment in intervals of duration P (here two values P = 10 3 , 10 4 will be considered) and 
we measure the correlation function between the first and the last times in the interval, C(0,P). This correlation is 
strictly related to the in-phase susceptibility, x' = (1 — C)/T, measured with an external magnetic field of frequency 
u> oc 1/P. In order to confirm this relation we show in the inset of Fig. [l] the susceptibilities measured in two reference 
experiments without any stop and with cooling and heating rates such that we perform P MCS at every temperature 
and then we change the temperature by AT = 0.02. The curves resemble very much the experimental ones. Compared 
to the AC measurements, simulations with the correlations have the advantage that one can reach higher values of 
P, i.e. lower frequencies. Note that, because of the precision required in this experiment, all the susceptibilities have 
been averaged over 100 samples of size 64 3 . In the inset of Fig. Ill some effects due to the finite cooling rates can also 
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be appreciated. In particular, one can see that the curves measured during the cooling processes stay a little bit above 
the corresponding curves measured during the heating process. This means that the system "accumulates" part of 
the relaxation work done in the low-temperature phase. The cooling and heating curves merge together only when 
the system comes back to the high-temperature phase. This phenomenon is present also in experimental spin glasses, 
even if with a much smaller intensitynJ. 

The interesting information can be obtained once we perform a long stop in the spin-glass phase. Here we stop the 
cooling at a temperature T* — 0.8 for t w — 100 • P MCS (i.e. 100 intervals) and measure x' from C(t w ,t w + P). The 
system relaxes and \'~ X re f becomes negative. However, when we start again cooling the system, the susceptibility does 
not recover completely the reference curve and it always remain with smaller values. Note that at the experimental 
levetm the susceptibility recovers the reference value very rapidly (around T ~ 0.65 -j- 0.7 on the scale of Fig. |l|) and 
that the apparent rapid increase just below T* in Fig. |l| is due to the very zoomed y axis and it has in fact a slope 
of order one. Moreover the curve followed by the data for T < T* does not seem to depend on P and if we consider 
the relative difference, (x' — Xref)/Xref> we would obtain that the convergence towards is still slower, due to the fact 
that both susceptibilities are decreasing with T. 

During the heating process the system stay on the same cooling curve and it does not show any strange effect near 
T*. For T > T* it finally recovers the reference curve and here is where we observe the largest dependence on P. 
The temperature where x' becomes comparable with x' le f strongly decreases with increasing P. Nevertheless for the 
times we have access to (P = 10 3 , 10 4 ) this temperature is larger than the critical one T c ~ 0.95 and in the limit of 
large times, P — ► oo, it can converge to both T c or T*. 

In conclusion we can assert that the three-dimensional Edwards- Anderson model does not show, on the time scales 
we have access to, the strong memory and chaos effects real spin glasses show. 

The same numerical experiment we have presented in this Section have been recently done also by Komori et al. (see 
Figure 7inref.M). In that fi gure it is shown that some kind of memory and chaotic effects are found when measuring 
the out-of-phase component of the AC susceptibility. But that figure and all the subsequent authors discussion based 
on it are inconclusive for the following reasons. They use a magnetic field that oscillates too fast (P — 160). As a 
consequence the effective critical temperature is very high, T° > 2.5 ~ 3T C . Moreover they use a cooling rate which 
is 10 times higher compared to experimental protocols. Each point in their figure is a measure over a single cycle of 
the field and a single dynamical history averaged over many different samples^]. Finally they claim to see in their 
figure 7 a merge of the susceptibility data to the reference curve, which is far from evident without any zoom of the 
interesting region. Note that both susceptibilities (x and Xref) goes to zero when T — > and then we also expect 
X ~ Xref to become zero. One should check that the relative difference is going to zero faster in order to claim for the 
presence of chaotic effects. In our study we had to increase the precision of more than 2 orders of magnitude (note the 
y-axis scale. in Fig. |l|) in order to discern the effect. The deceiving result is that, if the experiment corresponding to 
figure 7 oft 4 ] is done with slower cooling rates (unfortunately, such results were not shown intil) , one does not probably 
observe any trace of rejuvenation or memory due to the strong cooling rate effects. This is definitely different from 
what experiments show. 

B. Temperature cycling experiments 

In another set of very interesting experiments the temperature is changed according to the following scheduling: 
t w i seconds at T±, then t W 2 at Ti and finally t W 3 at T\ again. After that the TRM decay is measured. Depending on 
the sign of AT = T2 — T\ the system responds in different ways. For AT < an effective waiting time, , in the 
TRM decay can be defined and it is a monotonic function of AT such that i£f = t w \ + t W 2 + t w a for AT = and 
= twi+t W 3 for |AT| large. For AT > the TRM decay follow a more complicated law and it can not be described 
just by an effective waiting time. Nevertheless is always possible to define a correlation time, which turns out to be a 
monotonic function of AT, taking the same value as before for AT = 0, but converging to t^f = t W 3 for large |AT|. 

In numerical simulations the effective waiting time (or correlation time) can be estimated from the decay of the 
correlation function. In the inset of Fig. || we show the correlation measured in a cycling temperature experiment 
(with T\ = 0.7 and T2 = 0.9) where the relative times are similar to those used by experimentalists (t w i = 10 4 , 
t W 2 — 10 2 and t W 3 = 10 2 ). The reference curve with AT = will always refer to data measured at fixed temperature 
T\ and with a waiting time t w = t w \ + t w s. As it is clear from the data the temperature cycle does not affect at 
all the decay of the correlation function and the system does not seems to be re-initialized. This effect is still more 
drastic if we use a different scheduling in order to amplify it. In the main body of Fig. |^ we report the correlations 
measured with t w i = 5 • 10 5 , t W 2 = 10 6 and t w ^ = 5 ■ 10 5 , with T\ = 0.5 and T2 = 0.7 or 0.3. From the data it is clear 
that the effective waiting time is increased in both cases in contrast to what it is observed in experiments, being the 
time spent at the higher temperature T2 = 0.7 much more effective in terms of the thermalization process. 
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t 

FIG. 2. Correlation relaxation in temperature cycling experiments. In the inset we have used temperatures and time scale 
ratios similar to the experimental ones (see text for details) and we do not see any difference with respect to the reference 
curve. In the main part we show the results for a large perturbation both in temperature, AT = ±0.2, and in times, 
t W 2 = t w i + t w s — 10 6 . It is clear that for both positive and negative cycles the system is more thermalized with respect to the 
reference system where t w = 10 6 . 




We expect the effective waiting time to be more or less related to the size of thermalized regions in the system. 
Once we have characterized how this size changes under a temperature cycle, we can also study how the internal 
structure of these thermalized regions is modified by the cycle. In order to do this we exploit the generalization 
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of the fluctuation-dissipation relation to the off-equilibrium regime discussed in Section |Hl| In Fig. |3| we show the 
off-equilibrium susceptibility versus the correlation for the same experiments reported in Fig. || and described in the 
previous paragraph. The temperature cycle, even when it is very long, does not seem to affect the response of the 
system. In the quasi-equilibrium regime the slope of the curve x versus C gives the temperature in the thermalized 
regions. This temperature does not seem to change even if the system spends a lot of time in a different temperature 
T%. When it comes back to T% it rapidly seems to recover the configuration corresponding to temperature T\. 

C. Temperature shift experiments 

In order to better understand how the time the system spends at a temperature T\ can influence the thermalization 
process at a different temperature T2 we have performed a series of temperature shift experiments. Here the scheduling 
is the following. After t w \ MCS at temperature 2\ we set T — T2 and we immediately start measuring the correlation 
and the response to a small external field. In the present study T 2 = 0.5 and T\ — 0.7,0.5,0.3, the second case being 
considered as a reference curve. Moreover the waiting times, t w \ — 139, 10 3 , 10 5 , have been chosen such that is 
constant. It is known that in the Edwards- Anderson model the dynamical correlation length grows .as. a power law of 
time, £ oc t x ^ z ^ T \ where the dynamical exponent z(T) is inversely proportional to the temperatureE3c3. So our choice 
for the waiting times would correspond to thermalized regions of similar sizes. 




10 10 2 10 3 10 4 10 5 

t 

FIG. 4. The decays of the correlation function at temperature T = 0.5, after the thermalizations shown in the legend, have 
the same effective waiting time (as can be seen in the inset where we have rescaled the curves by means of simple multiplicative 
factors). The temperatures and the waiting times (see legend) have been chosen such that the size of thermalized regions, 
£(t, T) oc t aT (with a = 0.16 lEZJ and a prefactor of order 1), is comparable. 

In Fig. |^ we see that the effective waiting time generated by the three different scheduling is very similar. The 
three curves can be perfectly collapse in the aging regime by simply multiplying them by a constant (sec inset in 
Fig. This means that the effective waiting time is essentially given by the size of the thermalized regions, which 
has been chosen to be equal in the three experiments. The difference between the three curves in Fig. ^ comes from 
the quasi-equilibrium part C s t(t), which decays in a different way. 

The following natural question concerns the configuration of the system up to length scales of the order of 
Through the measure of the fluctuation-dissipation ratio we can estimate the effective temperature of the system 
on those length scales. In Fig. || we shown the results for the three experiments. The reference data, in the quasi- 
equilibrium regime, perfectly stay on the line (1 — C)/0.5 as it should. The other two data sets, because they have 
been thermalized at temperature T\ and then let to evolve at a different temperature T2 = 0.5, fall in between the line 
(1 — C)/Ti and (1 — C)/0.5, showing that the system temperature is changing from Tj to T 2 . However the interesting 
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point to note is that the change is very different in the two cases. For T\ — 0.3 (uppermost curve in Fig. |5|) the system 
responds with an effective temperature very similar to T2 = 0.5, given by the slope in the quasi-equilibrium regime. 
While for T\ — 0.7 (lowest curve in Fig. ||) the effective temperature is perfectly compatible with Tj = 0.7. In general, 
we would say that if an Edwards- Anderson model is thcrmalized to a temperature T\ and measurements are done at 
a different temperature T2 the response of the system will be dominated by the higher temperature. This is the only 
effect asymmetric in AT we,have found in all the numerical experiments performed with temperature changes. 

Recently Bernardi et alES have proposed the following scaling for the susceptibility x(tw,t w +t) = x(t;(t w ), 
where £(i) is the dynamical correlation length, defined above. Our data for the susceptibility which are measured on 
larger time and temperature scales than refE3, do not fit that scaling. 
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FIG. 5. FDT measurements taken at a temperature T = 0.5, after the thermalizations shown in the legend. They show that 
in the quasi-equilibrium regime (i.e. in the thermalized regions) the effective temperature of the system is always the maximum 
between the one used to thermalized the system and the one used to take measurements. 



V. MEMORY AND CHAOS IN AC RELAXATIONS 



In this Section we present a detailed investigation of memory, rejuvenation and cooling rate effects in the EA model 
doing AC susceptibility numerical experiments. In what follows we will use indistinctly the words chaotic, rejuvenation 
or reinitialization to indicate the presence of new relaxational processes which have been driven off-equilibrium when 
a perturbation is applied. There is a large list of AC experiments in the presence of external perturbations such as 
temperature field variations or magnetic field variationsBE3. The experimental setting is as follows. The system is 
quenched to a low temperature (ranging between 0.6 and 0.9 times the value of T g ) and the slow decay of the AC 
susceptibility recorded. After a time ti where ujt\ is such that the AC susceptibility has not totally decayed to its 
asymptotic value a perturbation is applied. To have an idea, the time t\ is such that both x' and x" are still between 
5% and 20%, of the whole decay, above their asymptotic large t value. This corresponds to typical values of ujt\ 
ranging from 100 to 2000. At this time there is a sudden perturbation (for instance, a change in temperature or field). 
After a time interval £2 which is of the same order as t\ the perturbation is switched off. In the present study and for 
sake of simplicity we have taken t2^t\. All the times t and frequencies lj we used in the numerical experiments are 
such that the scaling uit is satisfiedEd. 

In the presence of rejuvenation or chaotic effects one generally observes strong reinitialization of the AC susceptibil- 
ities corresponding to processes which have been driven off-equilibrium as consequence of the perturbation. Having in 
mind the previous experimental setting we have considered the following different types of perturbations: temperature 
changes, magnetic field changes and quenched disorder changes. 
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A. Temperature changes 



We quench the system to a temperature below T c and apply an AC magnetic field measuring the AC susceptibility. 
At a given time t\ we suddenly change the temperature to T + AT and still measure the AC susceptibility. Then, 
after a time interval equal to t\ we reset again the temperature to its original value T. In the presence of chaotic 
or rejuvenation effects we expect that a sudden change in temperature will reinitialize some relaxational processes. 
In Fig. ^ we show the results for L = 64 by measuring relaxation at T — 0.6 and making two jumps in temperature 
AT = ±0.1 at times t\ — 10000 and t% + = 2t\. The jump in temperature is then applied when a large part of 
the AC susceptibility is still relaxing like in the experimental setting. Note that for a positive temperature change 
AT the AC susceptibility stays above the reference curve at the temperature T + AT. This means that the effective 
waiting time after the positive jump is smaller than that of the reference curve at higher temperature. For a negative 
temperature jump —AT the AC susceptibility stays below the reference curve at temperature T — AT. This means 
that the effective waiting time has now increased and relaxation at T — AT has benefited from relaxation at the higher 
temperature T. The results of Fig. |6| show that the effective time i e ff during the interval of time ti when temperature 
has been changed is controlled by the same activated processes but with a different activation rate, 



(18) 



implying t e g > t 2 if AT > and vice versa. As comparison we also show a similar plot for t\ = 100000 when nearly 
all relaxation of AC susceptibility has taken place in Fig. 0. Again no traces of re-initialization effects in the AC 
susceptibility are observed. The main feature we can appreciate from the figure is that the positive jump and the 
time spent at T ± AT has increased the effective age of the system respect to the reference temperature T respect 
to the curve for the negative jump in agreement with Eq. (|l8|) . Note in the figure that, after the second jum p, the 
dashed line stays above-the continuous line. Note that this temperature dependence of the effective time (18|) is the 
one found in the REIVflo and here seems to behave quite well. 




log(t) 

FIG. 6. AC temperature variation experiment. The system is quenched at T — 0.6 and the AC susceptibility is recorded 
for field of period P — 100 and intensity ho = 0.1. At time ti = 10000 the temperature is changed by AT = ±0.1 and after 
the same time interval temperature is restored to its original value. As comparison we show the reference relaxation curves at 
temperatures T + AT = 0.7 and T — AT = 0.5 
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FIG. 7. The same experiment as in Fig. || but with ti = 100000. 
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FIG. 8. Plot above: The equivalent experiment of Fig. [] but with a larger size, a smaller field and a smaller frequency: 
L = 100, h — 0.01, P — 500. Plot below: The equivalent experiment of Fig. pi but with a 10 times smaller frequency P = 1000. 
Note that the perturbation in these two experiments and in Fig. M are applied when cot — 100. 



The picture which emerges from these figures is in agreement with all the results published up to now which point 
in the direction that there are no chaotic or rejuvenation effects in Edwards-Anderson spin glasses in the presence 
of temperature changes. One could argue that there are several factors which induce the absence of any trace of 
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rejuvenation. Among them: 1) The intensity of the field which is bigger than in experiments; 2) The smallness of 
the period of the oscillating field which covers at most nanoseconds when compared to real experiments (of order of 
seconds) and 3) The size of the system which is large enough. In Fig. || we show the same results as Fig. g for a field 
ho = 0.01, a larger size L — 100 and a larger period P = 500 such that uit takes the same value, so we are in the same 
time scale according to the results of ref.ES. Due to the smallness of the probing field the signal is now much more 
noisy so we show the in-phase susceptibility. The smallness of the absolute magnitude of the time scales involved in 
numerical experiments (reason 2) above) is usually advocated as the main source of discrepancy between numerics 
and experiments. The lower plot in Fig. g shows \' for the same size and field as Fig. ^ but with a frequency 10 times 
smaller. The conclusions are exactly the same: cooling rate effects are important and no trace of reinitialization after 
the temperature variation is observed. 

B. Field and coupling distribution variations 

To make evident how much this absence of chaos or rejuvenation is indeed an intrinsic effect to temperature changes 
we have done the same experiment with a different type of perturbation. Instead of changing the temperature, we have 
applied a perturbation which is well known to be .chaotic from equilibrium studies. Examples qf-such perturbations 
are: 1) a change in the uniform magnetic fieldOa and 2) a change in the couplings distribution!^. 

Concerning the first type of variation there have been several experiments which reveal how reinitialization occurs 
under a DC magnetic field changea. The experimental setting is the same as that shown previously but now the 
perturbation is to apply a DC magnetic field after t%. So the system is quenched at zero DC field and at time t\ the. 
DC field is switched on. After a time interval t% — t\ the DC field is set to zero again. In laboratory experimentst2l 
the intensity of the applied DC field must be larger than the amplitude of the AC field for the AC field to probe the 
response of the system after the DC perturbation. The intensity of the probing AC field is typically smaller than one 
Oersted and the intensity of the AC field much higher (between 5 and 10 Oersteds). So, typically the intensity of the 
DC field is 10 times or even more larger than the probing AC field. Nevertheless, in the numerical experiments we 
have a problem. The intensity of the AC field cannot be arbitrarily small, otherwise we have a too small and noisy 
signal. Consequently, if the perturbing DC field is chosen 10 times larger than the AC field then the resultant field 
(AC+DC) will be very large and drive the system out of the linear response regime. Moreover, non-linear effects will 
be much enhanced because of the non-linear coupling between the perturbation (the DC field) and the probing field 
(the AC field). A way to avoid this is to apply a perturbing DC field which does not couple with the probing field, 
for instance, a staggered DC field. Now the perturbation is given by a new term bJi in the Hamiltonian 

v 

5H = -h DC J2 e i (T i > ( 19 ) 
i=l 

where ti are quenched random variables which may take the values ±1. Another equivalent procedure would be to 
apply a uniform DC field as perturbation and measuring the AC susceptibility corresponding to the response to a 
staggered probing AC field. The results are shown in Fig. 0. Note that strong re-initialization is seen after perturbing 
the system in agreement with the known result that finite-dimensional spin glasses are chaotic against magnetic field 
changescfl. 

A similar result is found by considering the other type of perturbation. In that case we measure the AC susceptibility 
after quenching at temperature T. At t\ we take a percentage r (0 < r < 1) of the couplings and reverse its sign 
Jij = —Jij. After a new interval t\ we re-put the original couplings again. The results are shown in Fig. [l^ for 
t\ = 10000 10 = 0.01 and r = 0.05, 0.1 (corresponding to 5% and 10% of changes in the couplings respectively). Note 
the presence of strong and clear re-initialization effects in agreement with the fact that such a perturbation is chaotic. 

We may conclude this Section saying that while there is clear trace of chaotic behavior in the presence of field or 
couplings changes, there is absolutely no trace of re-initialization effects below the spin-glass transition in the case of 
temperature change experiments. This may be due to the presence of cooling rate effects in three dimensional Ising 
spin glasses stronger than those measured in laboratory experiments. 
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FIG. 9. DC Field variation experiment. The system L = 100 is quenched at T = 0.6 and the AC susceptibility is recorded 
for a probing field of period P — 100 and intensity ho — 0.1. At time ti = 10000, a DC field is applied for a time interval 
t2 — ti. After ti + t2 the field is switched off. The intensities of the DC fields are hoc—OA (triangles), 0.6 (circles) and 1.0 
(squares). The diamonds correspond to the reference curve without perturbation. 
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FIG. 10. Couplings variation experiment. The system L = 100 is quenched at T = 0.6 and the AC susceptibility is recorded 
for a probing field of period P — 100 and intensity ho — 0.1. At time ti = 10000, a percentage r of the couplings change sign 
and the AC susceptibility is recorded for a time interval t2 = ti. After ti + t2 the couplings take again their original values. 
The intensities of the perturbation are r = 0.05 (squares) and r = 0.1 (circles). The diamonds correspond to the reference 
curve without perturbation. 
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VI. OUTLOOK AND DISCUSSION 



Two years ago some experimentalists from the Saclay and the Uppsala groups measured on a spin glass sample 
very strong memory and chaos effectsa. Their results are really impressive and show unambiguously how important 
are these effects in real spin glasses. By investigating the Edwards- Anderson model in three dimensions we have tried 
to reproduce numerically their findings, but we have obtained results pointing in the opposite direction. Temperature 
variation experiments in spin glasses are nowadays one of the most puzzling results in the field. It is unclear which is 
the final theory, which may naturally account for these results. It is not easy to explain, from the point of view ftfrJrhf 
droplet modeE3, how reorganization of domains can account for the memory effects observed in the experiments&DLZl. 
Already in equilibrium theory, the question whether the spin-glass phase is chaotic against temperature changes 
is far fr om ibeing settled. Extensive numerical work does not show any clear evidence of chaos for temperature 
changesESE3c2l. 

Regarding off-equilibrium dynamics the situation is similar. Very precise numerical simulations by Rieger et al.E3 
show that correlations between equilibrium configurations at different temperatures are big and the corresponding 
overlap length grows algebraically in time without any tendency to saturation within the simulated range of times. The 
same.-epnclusions hold for TRM numerical experiments with temperature change protocols£3. Very recently, Komori 
et al.E-H have presented a detailed study of the two-time correlations in the presence of temperature change variations. 
These should be essentially equivalent to the present results because AC experiments probe the quasi-stationary aging 
regime where fluctuation-dissipation makes responses and correlations equivalent. Their conclusion is like ours: no 
rejuvenation effects are found and cooling rate effects are very strong. The only claimed evidence for rejuvenation and 
memory effects is the Figure 7 in that reference which is unconclusive as we have explained in Section fVL Contrarily, 
the results for DC field variation experiments resemble quite much those of experiments and, together with coupling 
distribution variations, show that chaotic effects manifest as re-initialization effects in the AC susceptibility. 

What is the origin of this discrepancy between experiments and simulations? We indicate three possible reasons: 

1. The size of the system is too small. 

2. The intensity of the probing field is too large. 

3. The absolute magnitude of the time scales in numerical experiments are too small. 

As we showed in Fig. || improvement in these limitations (larger size, smaller frequency, smaller AC field) does not 
seem to alter the final conclusions. Unfortunately we cannot still definitely conclude anything because probably our 
improvements (within the available computer capabilities) are too modest. Let us briefly comment about these three 
possibilities. Probably the finiteness of the size is the less important. A larger size diminislpapStaJjistical fluctuations 
but has no strong effect on the growing domains because they are known to be very smalBBEll Re garding the 
smallness of the field we have shown in Fig. || that a field 10 times smaller does not change the conclusions and for 
this smaller value of the probing AC field we are much closer in magnitude to the experimental setup. Moreover, 
field-cooling and zero- field coolirutesperiments, using the same values for the external DC field, show that the system 
is in the linear response regimeEilej. There are no deep reasons why things should drastically change for an AC 
probing field 100 times smaller to the one used in the present simulations. At most the intensity of the field can give 



The smallness of time scales involved (i.e. the fact that AC frequencies are too large and the waiting times too 
small) is the most serious reason and could definitively be the origin of the discrepancy. Short-time scales obviously 
imply short length scales of size £(P) or £(t w ) depending on the kind of experiment. If temperature changes would 
have no effect on small length scales but only on very large scales then chaotic effects could not be seen in standard 
numerical simulations. This statement corresponds to saying that the overlap length L(AT) for a typical temperature 
change AT is much larger than any probed domain length in the numerical experiment, but smaller than those probed, 
in the laboratory (where the chaos is clear). Let's try to quantify more this statement. In laboratory experiments^ 
typically length scales of the cpies of £ ~ 10 2 can be reached, while in numerical simulations we are restricted to 
length scales between 3 and 5e3cH. This discrepancy may or may not be a deep trouble depending on the value 
of the chaos exponent. Unfortunately there is not a numerical estimate for the chaos exponent (because .chaos in 
temperature has never been observed!!!) and we only have an estimate from domain- wall scaling argumentscj which 
give L(AT) = (AT)~^ with £ = d s /2 — 6 where d 8 is the fractal surface of droplets and fl-is the thermal exponent. 
For d = 3, the value of C, must be larger than 1 and a reasonable value seems to be £ ~ 1.5cZJ so chaos in temperature 
should not be too small after all. Note also that for a DC field change the chaos exponent is even smaller, the overlap 
length being given by L(h) ~ /i~3,o. So, for similar values for the prefactor--one would expect stronger chaos in 
temperature than in a DC field. On top of that, the analysis by Bray and MoorecS for the value of the chaos exponent 
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C, shows that is consequence of the balance between a contribution coming from the surface of the droplets L~ and 
a contribution from the activation energy necessary to revert a droplet L e . The surface contribution comes from the 
inhomogcncities in the couplings on the surface and, actually, the same type of analysis should fully carry through 
when analyzing chaos for coupling perturbations. But in this last case we found strong re-initialization effects (see 
Fig. |Io| ) which are absent when changing the temperature. How this difference which we observe can be explained 
in the framework of the droplet model remains mysterious unless one advocates different prefactors for both overlap 
lengths (much bigger for temperature changes than for couplings changes) or that something special occurs in three 
dimensions^. With the present estimations for the £ exponent and if the prefactor for the overlap length is very large, 
chaotic effects, which experimentally appear on temperature changes such that AT ~ 0.1 T Cl may need in present 
simulations temperature changes of the same order of the absolute temperature. If this is the case then we have to 
wait for the next computer generation or to study the EA model in a situation such that larger length scales could 
be reached (e.g. in 4d or in 3d with next-nearest neighbours interactions). 

On the other hand our results imply that if the laboratory experiments were done at frequencies of 10 9 Hertz 
(instead of the typical 1 Hz measurements) then cooling rate effects would be restored and rejuvenation or chaotic 
effects disappear. Unfortunately, there are no experiments in the range 10 — 10 9 Hertz. Note also that, according to 
the general tut scaling, for these frequencies one should do measurements at very short times such-that cut is not much 
larger than 10 4 and the AC susceptibility has not completely relaxed. Still, the results of ref.ej, the experimental 
results ofl and all the numerical published data up to now show that the scaling tut (the equivalent of the scaling t/t w 
in two-time experiments) works reasonably well (with some slight deviations). If the scaling t/t w means something (as 
most of the present theoretical work suggests) then it is difficult to understand why no trace of re-initialization effects 
is observed in the smallest coarsening domains when the temperature is changed. Actually these re-initialization 
effects are found for DC magnetic field and coupling distribution changes. If some dynamical effects are completely 
absent for the small coarsening sizes this means that no numerical simulation in the last ten years has actually reached 
the asymptotic regime where connection to real experiments is possible and we are certainly missing somethingEZI. 

Concerning theory, it is difficult-to give a complete description on these effects in terms of compact excitations 
as proposed in the droplet modeled. Usually the asymmetric response of the spin glass against the sign of the 
temperature perturbation is explained in terms of an asymmetric overlap length L(AT). The problem is that in 
simulations L(AT) is apparently extremely large when compared to experiments and also symmetric, a question 
which is difficult to explain again if one does not appeal to the smallness of the time scales involved. If the necessary 
time scales to see the asymmetric effects in L(AT) exceed the experimental time scalesE?J then it is unclear how to 
reconcile any two among the three: theory, experiments and simulations. p, 

The explanation of memory and chaotic effects was originally explained in terms of a hierarchical picturetj. Actually, 
these effects are explained and reproducible in theJjjEEMEHfE 1 ], a model with several critical temperatures. The same 
phenomena is absent in the random-energy modelEH with a single critical temperature. Unfortunately, the GREM 
is a model without microscopic and spatial description so the connection with real spin glasses remains speculative. 
As we said previously the Jesuits we find here, and in particular the cooling rate effects, are very similar to those of 
the random energy modelEa. It is well known that the spin-glass phase in this model can be described in terms of 
a one-step solution. And what we see-ia our simulations is what is expected for a mean- field model with a one-step 
of replica symmetry breaking solutionEa. This solution is known to. describe the physics behind structural glassesa. 
Actually temperature variation experiments on structural glassesE3~E3 resemble our simulation results much more 
than what experiments on spin glasses do. Again one could claim that the one-step character of the effects we observe 
in our simulations are due to the smallness of time scales and if one increases the time scales (let us say by 6 orders 
of magnitude) much different results will come out. 

But there could be another and more natural explanation. Is the Edwards- Anderson model really a good spin glass? 
Is it possible that the Edwards- Anderson model has strong cooling rate effects and no rejuvenation effects at all? If 
rejuvenation or chaotic effects are present in the dynamics for a given perturbation this could be related to the fact that 
equilibrium properties are chaotic when such a perturbation is switched on. Analytical calculations in the mean-field 
version of the Edwards- Anderson .model (the Sherrington-Kirkpatrick model) confirming chaoticity for temperature 
changes are still inconclusiveotJ. A completely different scenario holds for field and coupling perturbations in 
agreement with the present simulations. Let us note also that real spin glasses are site disordered systems, a type p£ 
disorder not included in the Edwards- Anderson model. Recent experiments on the Kagome antiferromagnet latticeEZI 
reveal that there are not strong re-initialization effects after changing the temperature. Because that system does 
not include disorder at all one may wonder whether the influence of site disorder can. be important. Site-disordered 
spin glasses were not much considered in the past (for some works and references seeu) because they were thought to 
be less relevant to experiments than bond-disordered models but probably this is not true and a site disorder effect 
should be taken into account to explain experiments. To finish this collection of possible ways out to this puzzle let 
us point out also the possible role of the continuous character of the spins in a real spin glass as well as the effect of 
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chiralityEl This last effect and its importance in the description of the dynamics of spin glasses will surely see further 
developments in the forthcoming years. 

It is clear that we are facing a very difficult problem. If the smallness of time scales of the simulation is the 
final explanation for everything then a sensible theory for spin glasses must explain why temperature changes are 
fundamentally so peculiar when compared to other type of perturbations, i.e. why the prefactor for the overlap length 
L(AT) is so large. In this respect, experiments at much larger frequencies are necessary. At least, to see if cooling 
rate effects gradually change and re-initialization effects, in the presence of temperature variations, do systematically 
weaken. On the other hand, if the Edwards-Anderson model fails to explain a crucial result found in experiments 
then we must discover what ingredient is lacking in the original model and what are the consequences for our present 
knowledge of spin glass theory. In this direction, finding a microscopic model with spatial structure which presents a 
spin-glass phase transition and unambiguously shows memory, chaotic as well as absence of cooling rate effects would 
be welcome. 
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